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Abstract 

We present a novel cosmological model in which scalar field matter in a biaxial Bianchi IX geometry leads 
to a non-singular 'pancaking' solution: the hypersurface volume goes to zero instantaneously at the 'Big 
Bang', but all physical quantities, such as curvature invariants and the matter energy density remain finite, 
and continue smoothly through the Big Bang. We demonstrate that there exist geodesies extending through 
the Big Bang, but that there are also incomplete geodesies that spiral infinitely around a topologically 
closed spatial dimension at the Big Bang, rendering it, at worst, a quasi-regular singularity. The model 
is thus reminiscent of the Taub-NUT vacuum solution in that it has biaxial Bianchi IX geometry and its 
evolution exhibits a dimensionality reduction at a quasi-regular singularity; the two models are, however, 
rather different, as we will show in a future work. Here we concentrate on the cosmological implications of 
our model and show how the scalar field drives both isotropisation and inflation, thus raising the question 
of whether structure on the largest scales was laid down at a time when the universe was still oblate (as also 
suggested by lIllEll^l)- We also discuss the stability of our model to small perturbations around biaxiality 
and draw an analogy with cosmological perturbations. We conclude by presenting a separate, bouncing 
solution, which generalises the known bouncing solution in closed FRW universes. 
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I. INTRODUCTION 



At the core of most of theoretical cosmology lie the assumptions of homogeneity and isotropy. 
These were originally motivated by the cosmological principle and the mathematical tractability 
of the resulting FRW models. However, the observed universe is obviously neither homogeneous 
nor isotropic, so these symmetries can only ever be approximate. Thus the question arises as 
to which assumptions we can relax whilst maintaining analytical tractability. We choose here to 
consider more general, anisotropic cosmologies, whilst keeping the assumption of (large-scale) 
spatial homogeneity. The resulting class of cosmologies is collectively known as the Bianchi 
models. 

Scalar fields are ubiquitous in theories of high energy physics. The Standard Model of particle 
physics postulates the scalar Higgs particle, and Superstring Theory and string-inspired models 
motivate a plethora of moduli fields from compactifications, dilatons, radions etc. Scalar fields are 
also of interest in cosmology, as they can drive periods of accelerated expansion of the universe 
in a relatively straightforward and plausible manner. Cosmological scalar fields could therefore 
account for the posited period of inflation and the apparent present period of A-domination. 

It is therefore natural to consider the effects of a scalar field dominating the dynamics of a 
Bianchi model. However, the phenomenological success of FRW models suggests that we should 
only consider the Bianchi classes that allow FRW universes as special cases. This narrows the 
phenomenologically interesting Bianchi types down to I, V, VII^ and IX. From previous work 
m [51 we are particularly interested in closed models, which are still consistent with observations 
||6l. The only closed Bianchi type that also allows for FRW subclasses is Bianchi IX, so we shall 
consider the dynamics of a scalar field in a Bianchi IX model [H [8l HI. Bianchi IX models are 
known to have very complicated dynamics, exhibiting oscillatory singularities and chaos (mix- 
master) [ni[l2l[l3l[14l[T5l[l6l[l7l. In this light, it is even more intriguing that the biaxial 
model we will consider is so well-behaved. 

It is interesting that our solution provides another example of a regular solution with axial 
symmetry; for others, see Pitrou et al [[B 13 and Senovilla and collaborators [fTSl [T9l l20l |2T| . 
In the different context of inhomogeneous but axially symmetric cosmologies, iflSl [T9l l20l |2T]| 
also consider the important concept of geodesic completeness of a spacetime. Completeness is 
a criterion for the physical significance of a solution that is more stringent than mere geometric 
regularity, which we will also need to address later. 
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Often inflation is invoked to justify the assumptions of (acausal) homogeneity and isotropy. 
However, this then raises the question of how sensitive inflation itself is to the rather unnatural 
initial conditions of a homogeneous initial state. For these considerations, the interested reader 
may refer to, for instance, [|22|. Here, however, we consider a Bianchi geometry as our starting 
assumption, without invoking inflation in order to justify it. For recent related literature concerning 
the case where a vector field is present to drive a phase of anisotropic inflation, see for example 
1231121. 

This paper is organised as follows. We begin with a brief introduction to Bianchi models 
in Section [11} In Section [111] we construct a generic Bianchi IX model dominated by scalar field 
matter and then specialise to the biaxial case in Section |lVj (An alternative form of the Einstein 
field equations for this system, using the 3 + 1 covariant approach, is presented in the appendix.) 
We demonstrate scaling behaviour of solutions, and present a particular series solution, in which 
one radius goes to zero at the Big Bang (pancaking). We also consider geodesies through the Big 
Bang. We then present a realistic cosmology based on our model in Section [V] and, in Section |VI[ 
the stability of such a model about biaxiality. We finally present a separate, bouncing solution in 
Section [VIll before we conclude in Section [Villi 

II. BIANCHI UNIVERSES 

Bianchi universes are spatially homogeneous and therefore have a 3-dimensional group of 
isometrics G3 acting simply transitively on spacelike hypersurfaces. The standard classification 
follows Bianchi's (1897) classification of 3-parameter Lie groups ll25l . 

We adopt the metric convention (H ). Roman letters a,b,c... from the beginning of the 

alphabet denote Lie algebra indices. Greek letters /i,v,(7... label spacetime indices, whereas 
Roman letters i.j.k... from the middle of the alphabet label purely spatial ones. 

The isometry group of a manifold is isomorphic to some Lie group G and the Killing vectors 
obey ^v] = C^yi^^ where the C^y are the structure constants of G, so Lie groups can be used 
to describe symmetries, in particular, isometrics. 

When studying symmetries, an invariant basis is often useful. This is a set of vector fields 
each of which is invariant under G, i.e. has vanishing Lie derivative with respect to all the Killing 
vectors such that 

[^M'^v]=0. (1) 
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Such a basis can be constructed simply by imposing this relation at a point for some chosen set of 
independent vector fields and using the Killing vectors to drag them out across the manifold. The 
integrability condition for this set of first-order differential equations amounts to demanding that 
the C^y be the structure constants of some group. The invariant vector fields in fact satisfy 

[Z^,Zv] = — C^yX(j. (2) 

Denoting the duals of the by (0^, the corresponding curl relations for the dual basis are 

dco^ = ^C^cjtC0'' A(o\ (3) 

Because the are invariant vectors, the metric can now be expressed as 

ds^ = g^yCO^Oi^, (4) 

for some 

Bianchi models can be constructed in various different ways. They are based on properties of 
a tetrad {e,} that commutes with the basis of Killing vectors {(^/} which generate the symmetry 
group 

[eoy=0. (5) 

Here we have made use of the fact that for homogeneous cosmologies there exists a preferred 
foliation of the 4-dimensional spacetime into a product spacetime with a timelike Killing vector 
and a 3-dimensional group of isometrics G3 acting simply transitively on spacelike hypersurfaces 
such that the timelike Killing vector commutes with the spacelike Killing vectors. Therefore we 
change from spacetime indices to spatial indices. We also adopt the comoving gauge whereby 
the timelike basis vector is taken as parallel to the unit normal to the surfaces of homogeneity. 
In general, the structure constants of the symmetry groups of these different surfaces of constant 
time are time-dependent [fTSl l26l l27ll . However, from ([5]), it follows that between homogeneous 
slices, the structure constants are preserved up to a time-dependent linear transformation. One is 
free to shift the time dependence between the structure constants and the metric components by 
adjusting the time evolution of the tetrad. We choose to put all the time variation in the spatial 
metric components such that 

ds^ = dt^ - Yki{t){e'^{x)dx'){e'j{x)dx-'), (6) 
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where ef (x) are the one-forms inverse to the spatial tetrad which have the same structure constants 
C^jj as the isometry group and commute with the timelike unit normal eo to the surfaces of ho- 
mogeneity: eo = dt and e,- = ejdj, such that [ei.ej] = C^j^k and [eo,e/] = 0. This is the approach 
we will adopt for constructing the metric later, as the Einstein field equations become ordinary 
differential equations in the metric components yij{t). (Other approaches are based on the auto- 
morphism group of the symmetry group [|28l |29l or put the time-dependence in the commutation 
functions of the basis vectors [[30ll .) 

We will now briefly review the Bianchi classification of G3 group types. The spatial part of the 
structure constants C^y can be decomposed into irreducible parts as follows 

C]j = £ijin'^ + at 5f - a J df (7) 

for n'^ symmetric. The Jacobi identities can then be rewritten as 

n'^aj = 0. (8) 

Without loss of generality we can choose the tetrad so as to diagonalise «,j = diag (ni,n2,n3) 
and to set a' = (a, 0,0) which reduces the Jacobi identities to nia = 0. This then allows one to 
classify the possible Lie groups. We can define two broad classes of structure constants according 
to whether a = (Class A) or not (Class B). See Table|I]for more details concerning the individual 
Bianchi types specified by the different options for the n^j matrix entries. 

If one were to put the time dependence in the structure constants instead, one would of course 
have to show that the classification type must be preserved by the evolution equations for n{t) 
and a{t). It turns out that it is a generic property of the Einstein field equations that they pre- 
serve symmetries in initial data within its Cauchy development, so that the classification is in fact 
independent of which approach one chooses OTll . 



III. TRIAXIAL BIANCHI IX MODEL 



We now consider Einstein-Hilbert gravity in a generic Bianchi IX model with a minimally 
coupled scalar field. For generality we also include a cosmological constant term. We therefore 
start with the action 



1 

2k 



1 



i? + 2A)--V^(^)V^(^) + \/(<^) 



(9) 
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TABLE I: Bianchi model classification. Those containing FRW models as special cases are: Bianchi IX 
(closed); Bianchi I and Bianchi VIIq (flat); Bianchi V and Bianchi Vllh (open). 

where R is the Ricci scalar, ^ the scalar field and V{^) its potential, which we will for simplicity 
assume to be that of a simple massive scalar field, V(0) = jnp-(j)^. Variation of the action with 
respect to the scalar field yields conservation of energy-momentum 

Vjur^'' = 0, (10) 

for the usual scalar field energy-momentum tensor 

T^tv = M;v-8nv (i0'^0;p -y (</))) . (11) 

Variation with respect to the metric yields the Einstein field equations 

G^y=KT^y+Ag^y. (12) 

We are working in Planck units where c = h= I and we will set K = StiG to unity eventually. 
We choose to express the Bianchi IX metric in the form 

ds^ ^dt^-^Yij{t)(o'coj (13) 
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for yij{t) = diag {R\{t),R\{t),R\{t)) and invariant 1 -forms co' given by 

(0^ = dx + sinydz, 

CO = cosxdy — s'mxcosydz: 

(0^ = smxdy + cosxcosy dz, 



(14) 



with corresponding Killing vectors (these just correspond to the rotations of the 3-sphere, i.e. the 
Clifford translations) 



and invariant basis 



^1 

^3 

^1 

X3 



secy COS z (9^ + sinz (9^, — tanj COS z (9^, 
— sec J sin z + cos z dy + tan y sin z (9^ , 



(15) 



sinxtan3;(9x-fcosx(9y — sinxsec3;(9j;, 
— cosxtany (9x + sin.x;(9y + cos.x;sec);(92 



(16) 



taken from [|32l . (We will consider an alternative set of Killing vectors derived within the Con- 
formal Geometric Algebra framework [|33l |34l related more obviously to an {x,y,z) coordinate 
system in a future work.) 

Using the above expressions for the generating one-forms and expanding out we get the fol- 
lowing non-zero metric components for a triaxial Bianchi IX universe (where we have dropped 
the explicit time dependence of the scale factors Ri{t) for the sake of brevity): 



gtt 



Syy 



1 n2 



r(i?2COS^A:-fi?3 sin^x) 



Szz = sin^j-f- (i?2sin^jc-f-i?3Cos^x)cos^y], 

gxz = -5^? sin J, 

gyz = — 5(i?3 — i?2) sinxcosxcosj. 



Introducing the usual definitions for the Hubble parameters 



R, 



(17) 



(18) 



for the three different directions, we find that for the above metric the conservation of energy- 
momentum ( [To] ) corresponds to the Klein-Gordon-type equation 

m^(j) + {Hi+H2+H3)(j) + (p = 0. (19) 

The ^/-components of the Einstein field equations (for spatial index /) give three dynamical equa- 
tions for the Hubble parameters Hi, the first of which reads 

2Hi = H2H3 - 2Hj - H1H3 - H1H2 + A-KP 
5Rl 3Rl 3Rl 6 2 2 



nlnl^ nlnl^ nlnl nl^ nl' ^^^^ 

K2i\2 K^K^ /\j/\2 /Vj K2 

and the other two equations are obtained simply by swapping indices 1^2 and 1^3. The 
??-component of the Einstein field equations yields one Friedmann-type constraint equation 

-H1H2- H2H3 - H1H3 +A+Kp 

i?? i?T i?Q 2 2 2 „ 

^ p2n2 ^ p2p2 ^ p2p2 p2 p2 p2 ^ ^ 

K2i\2 /Vj/Vj K^K2 K2 

where p = + V {(j)) and p = — V {(j)) are the energy density and pressure of the scalar field 
matter respectively. Note that there are no spatial gradients of ^ in the expressions for p and p, by 
spatial homogeneity. The case of vanishing potential corresponds to a stiff fluid, p = p, which has 
been investigated in [[35ll36ll . amongst others. 

For different applications, it might be useful to recast these equations in terms of the averaged 
scale factor, i.e. the volume expansion, 

R = {RiR2R3)K (22) 



its associated Hubble factor. 



3H = Hi+H2 + H3, (23) 



and the shear Oij. For example, this would be useful for separating the contributions from the 
curvature and the shear. We display this form of the Einstein field equations in the appendix. In 
the case of the biaxial Bianchi IX model that we will consider shortly, however, we believe the 
parametrisation in terms of the different radii is clearer. When considering the cosmology at late 
times, however, we will find that our model has isotropised sufficiently such that we can describe 
our model by an effective FRW-model with a scale factor given by the volume expansion of the 
Bianchi model. 
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It has been known for a long time that the Bianchi IX model (or indeed any Bianchi model 
for which C\j = 0) is geodesically incomplete and exhibits a curvature singularity for perfect 
fluid matter (see, e.g. [|37ll ). The singularity occurs precisely when the volume of an invariant 
hypersurface goes to zero. The perfect fluid energy density can be shown to diverge as py^/ ~ y 
and the curvature invariants are singular as well, as ? ^ 0. However, it is the oscillatory character 
of the singularity that makes perfect fluid matter in a Bianchi IX universe model unsatisfactory as 
a cosmological model: the Big Bang is an essential singularity. Oscillations in the ratios of the 
different scale factors in ln{t) effectively show that there is an infinite history upon approaching 
the Big Bang, making it impossible to trace back to it. Conversely, coming out of an essential 
singularity to reach the observed universe is ill-defined for the same reasons. This also holds, in 
particular, for the biaxial case by virtue of the general theorem. We will see below that for a biaxial 
Bianchi IX model with scalar field matter, however, the picture will be qualitatively different. 

Scalar field matter is still unsatisfactory in the full triaxial Bianchi IX model, however, since 
it exhibits the same behaviour outlined above. As we show below, it is only when we impose 
axial symmetry, that a solution is possible in which all physical quantities such as energy density 
and curvature invariants remain finite at the Big Bang, and the universe extends smoothly across 
what is no longer an essential singularity into a well-behaved pre-Big Bang phase (though parity- 
inverted!). This then raises the question as to whether the axially symmetric case is stable to small 



pertubations in biaxiality, which we will address in Section VI 



IV. BIAXIAL BIANCHI IX MODEL 

We therefore now specialise to the case where two of the axes are equal, Riit) — R3{t), which 
leads to the simplified metric (again with the ^-dependence of the scale factors suppressed) 

/ 1 \ 

-\rI -iRjsiny 



-\rI 



(24) 



\0 -\R\%my -\{R\%m^y + Rlco^^y) ) 

Now there are only two dynamical equations, for the two non-degenerate Hubble parameters, 
and one Friedmann constraint. The first dynamical equation is 

IH2 + 3//| + - A = (^3^ - 4^ (25) 



and the other reduces to 

2Hx + 2Hl - Hi + 2H1H2 + Jc;? - A = - (^5^ - 4^ . (26) 

In the isotropic limit, these equations reduce to appropriate combinations of the usual Friedmann 
and acceleration equations, as required. The Einstein field equations further yield the Friedmann 
constraint 

Hi + 2H1H2 - Kp - A = -1 (^^ - 4^ . (27) 

This is straightforwardly seen to reduce to the standard Friedmann equation as 7?i ^ i?2 = ^3- 
The equations of motion for the matter content follow straightforwardly from energy-momentum 
conservation as above. In fact they are easily deduced from the triaxial case and read 

n?(^ + {Hi + 2H2) ^ + (f>=0. (28) 



We prefer the viewpoint from which the equations (25 )-(26 1 are dynamical equations for the Hub- 
ble parameters and regard the scale factors Rt as derived quantities, but these equations are obvi- 
ously equivalent to second order equations in terms of the radii. 

The simplicity of these equations suggests that a relatively simple solution should be possible. 
We are particularly interested here in solutions with definite parity, and defer discussion of solu- 
tions with indefinite parity to a future work. We show below that both odd-parity and even-parity 
series expansions exist around the 'Big Bang', which are therefore valid starting points for numer- 
ical integration, but first we briefly discuss the generation of a family of solutions from a given 
solution by scaling. 

A. Families of solutions related by scaling 



Given a solution to the equations (|25|)-([28]), a family of solutions is generated by scaling with a 
constant o and defining 

1 - - - T 

R^{t) = -Ri{ot), Hi{t) = oHi{ot), ^{t) = (l){ot), m = om, A = o^A. (29) 

(This scaling is in fact analogous to the one found previously in dU.) This scaling property is 
valuable for numerical work, as a range of situations can be covered by a single numerical integra- 
tion. Furthermore, many physically interesting quantities turn out to be invariant under changes 
in scale. This scaling property does not, however, survive quantisation, so one would have to be 
careful when considering vacuum fluctuations. 
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B. Odd-parity series expansion around the 'Big Bang' 



It is well known [[I5l[371 that pancake singularities (which we will identify with the 'Big Bang', 
and choose to happen at time t = 0), where one radius goes to zero and the other two remain finite 
as ? ^ 0, can occur in Bianchi models. This is already in some sense an improvement over the 
FRW case, where the singularity is pointlike. However, pancake singularities are known to occur 
in Bianchi I (as well as cigar singularities), whereas Bianchi IX is commonly thought generically 
to exhibit oscillatory singularities. We will now show that, with scalar field matter, the pancake 
singularity in the Bianchi IX model is not in fact a curvature singularity at all, because all physical 
quantities remain finite as ? ^ and extend smoothly into a parity-inverted universe for t <0. 

We assume that it is the non-degenerate radius Ri that tends to zero as ? — > 0. Close to the Big 
Bang, the linear term becomes dominant for Ri{t), whereas the other radii ^2(0 = ^3(0 t^^d to 
a non-zero constant, as does the scalar field ^{t). Indeed, for a series solution ansatz in which we 
assume oddness for Ri and evenness for the other functions, i.e. 

Ri{t) =t [ao + a2t^ + a4t'^ + ...) 
R2{t)=R3{t)=bo + b2t^ + b4t^ + ... , 

0(O=/o + /2?'+/4?V... 

(30) 

then the three dynamical equations ( [25] ), ( [26] ) and ( [28] ) allow us to fix the three series term-by- 
term, given the initial values gq = Ri (0), bo = RjiO) and /o = (^(0). The fact that this also satisfies 
the Friedmann energy constraint ( [27] ) then proves that this odd-parity series solution is a valid 
expansion around the Big Bang, which we can use as a starting point for numerical integration. 

Intriguingly, it turns out that the spacetime is non-singular insofar as the Riemann tensor is 
well-behaved at the Big Bang and so are all curvature invariants. In fact only three components of 
the Riemann tensor are non-zero at ? = 0: 

Rtyty = ^KmVo^O + l^bl - i, Rtztz = RtytyCOS^y, Ryzyz = -\blcOS^y. (31) 

The fact that there is no curvature singularity means that these series solutions can be contin- 
ued through to negative values of t. That is, R2 = R3 and are even for t <0, whereas Ri is odd 
i.e. negative for t < 0. This amounts to a parity inversion at the Big Bang. Instantaneously, as 
the hypersurface volume goes to zero at the Big Bang, the spatial hypersurface is 2-dimensional. 
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The fact that the universe can be momentarily 'dimensionally reduced' is interesting. Nonethe- 
less, even though the hypersurface volume goes to zero instantaneously, information about the 
evolution is encoded in the derivatives. If one adopts the string-inspired viewpoint [|38l |39l |40ll 
that the universe could be described by a 3-brane, the instantaneous conversion of 3-branes into 
2-branes seems problematic in type IIA and IIB String theory, but could be interesting from a type 
IIA/IIB String theory duality or AdS/CFT point of view. From a more conservative perspective, 
this process resembles Taub-NUT space ll37l l39l , which is a biaxial Bianchi IX vacuum solution 
that can be represented as a disc that evolves into an ellipsoid and back into a disc. It therefore 
also shows the same feature of dimensional reduction as our model, and moreover also does not 
have a geometric singularity during this collapse. In particular, it is thought to evolve from time- 
like open sections in a NUT region, via lightlike sections (called Misner bridges), to spacelike 
closed sections in the Taub region, back into timelike open sections in the other NUT region. This 
open-to-closed-to-open transition is not mathematically singular, but it is incomplete, as geodesies 
spiral infinitely many times around the topologically closed spatial dimension as they approach the 
boundary [|4ni42]| . This type of singularity is called 'quasiregular' in the Ellis and Schmidt classi- 
fication [[43]| (these include the well-known 'conical' singularities [|44l), as opposed to a (scalar or 
non-scalar) curvature singularity. For completeness, a possible parametrisation for the Taub-NUT 
metric is 

ds^ = Idtco^ - iRlico^)^ - iRj [(«2)2 + (©3)2] (32) 

as compared with 

ds^ = dt^ - \R\{o)^f - \rI [{03^f + {(o^f] (33) 

in our model. The interested reader may refer to ll37l . for example, for a discussion of Taub-NUT. 

It is worth noting that in deriving the above series solution, one can instead start with a general 
Taylor series expansion in the full triaxial case. Expanding around the Big Bang rather than any 
other point amounts to demanding that the constant term in the Taylor expansion for R\ be zero (in 
fact, this could be any radius, as Ri is not special in the triaxial case), whereas everything else is 
in principle undetermined. Imposing the dynamical equations, however, forces the two non-zero 
radii to be the same term-by-term for the expansion to be valid, i.e. a Big Bang-like expansion 
only works for the axisymmetric model; there is an essential singularity unless we consider the 
special, biaxial case. The validity of the equations of motion further implies the above mentioned 
even form for R2 = R3 and ^, and odd behaviour for i?i . 
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Since the existence of the above series expansion around the Big Bang demonstrates that an 
axisymmetric pancake singularity is a valid starting point for numerical integration, we therefore 
solved the biaxial equations numerically subject to the appropriate boundary conditions. We found 
good agreement with the series solutions within the range of their validity. The large parameter 
space of this model admits both viable (in the sense of compatible with present observations) 
and unrealistic (in the sense of incompatible with present observations) cosmologies. We defer 
seeking realistic model parameters and displaying the numerical results until Section|V]and instead 
concentrate next on addressing the issue of geodesic completeness. 



C. Behaviour of geodesies 

As we have shown, our model has no curvature singularity at the Big Bang, at which all physical 
quantities remain finite. Thus, at the level of the evolution equations, the model is well-behaved. 
We now consider the question of geodesic completeness of our model to determine whether it 
possesses a non-curvature singularity at ? = 0. 

The geodesic equations are most easily obtained using the Lagrangian formalism, in which 

L=lg^yx^x^ (34) 

is varied with respect to the coordinates [x^] = {t,x,y,z)', here a dot denotes a derivative with 
respect to some affine parameter A and the factor ^ is included for later convenience. Inserting the 



biaxial metric ( 24 ) into L yields 

L = 4/2 - Rjx^ - 2R\ sinjiz - Rlf - {r\ sin^y cos^y) z^. (35) 

Since L is independent of the x- and z-coordinates, the corresponding Euler-Lagrange equations 
yield two conserved quantities and according to the relations 

-lR\{x + ^myz) = K^, (36) 
-2[Rlsmy{x + smyz)+Rlcos^y-z] = K^, (37) 

which may be solved for x and z to yield 

KzR\^m.y-Kx{R\^\TL-y + R\co^^y) 



IRjRlcos^y 

2i?2 cos^y 
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(38) 
(39) 



Substituting these expressions back into the Lagrangian, the Euler-Lagrange equation for y then 
reads 

ARI cos^ y (^2y + 2 = [K^ + siny - 2K^K, + K^K^ cos^ y. (40) 

Let us first consider null geodesies for which g^j^yX^x^ = = L/4. To begin with, we wish only 
to show that some geodesies extend smoothly through the Big Bang, and therefore we can choose 



a convenient form. Inspection of equations ( [38] ) and ( [39] ) shows that we cannot set either of i and z 
alone to zero by choosing suitable constants, but we can consider the case where both vanish (this 
requires the choice Kx = K^ = and essentially corresponds to motion in the y-direction). In this 
case the Euler-Lagrange equation for y simplifies to d {R\y) /dX =0 which can be immediately 
integrated to yield 

Rly = Ky, (41) 

where Ky is another constant of the motion. The constraint that the geodesic is null now reads 
— R^-)!^ = 0, which gives 

^ = ^ = ±-. (42) 
dt t R2 

This is reminiscent of the Friedmann case, but one must now bear in mind that in our model R2 is 
finite and non-zero at the Big Bang in contrast to the scale factor in an FRW model. Thus, there 
is no longer a difficulty in integrating through the Big Bang; rather, it is trivial. Indeed, looking at 
the Euler-Lagrange equation for y directly in this case, one has 

i?2j- + 2^/3) = 0, (43) 

which amounts to j' = near the Big Bang where R2 has a minimum, also suggesting that such a 
light ray goes straight through the Big Bang. Thus photon motion in y is effectively undisturbed 
by the pancaking. 



Let us now consider a photon moving in the x-direction. It is easiest to see from (40) that if a 



photon starts at y = it will stay there. Hence for = and y = 0, the results (38 ), (39 1 and (40) 
yield i: = 0, j = and 

(44) 

2R\ 

In this case, the constraint that the geodesic is null yields 4/^ = R\x^. Hence, close to ? = 0, where 
Ri is approximately linear in t, once has 

— = - = ±— oc ±-. (45) 
dt t Ri t 
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Integrating we obtain 

x{t) = ci\n{-t)+C2 or x{t) = C3\n{t) + C4, (46) 

for some constants c,, z = 1 . . . 4, which shows that at ? < we must have the former solution and 
at f > the latter. 

Let us illustrate these solutions by considering the specific case of a photon travelling in the 
positive jc-direction as time increases, i.e. j§ > 0. For ? < we have ^ = 7- and for this to 



be positive for negative t we must have ci < 0. This also yields Ri ——t from (45) which is 

C\ 1 



negative for ? < as expected. For f > (46) yields ^ = ^ > so C3 > and i?i ^ ^t. So 



it seems that the photon x-coordinate goes to +0° as ? ^ from below, reappearing just after 
t = dX X = —00. This is less problematic than it seems: the spatial sections in Bianchi IX are 
topologically 5^, so it is natural to think of the coordinates as angles that should be periodically 
identified. Furthermore, at ? = and therefore 7?i = 0, (|6]) yields vanishing proper distance between 
different values of x, which presumably means that this direction has shrunk to a point. Since our 
coordinates are angles, the result of the photon going to 00 appears to signify that the geodesies are 
infinitely spiralling (which is reminiscent of Taub-NUT) around a closed spatial dimension as it 
is collapsing. The proper distance traversed and time taken both, however, become infinitesimally 
small, so the whole process might still be finite. Thus it is not clear that anything singular has 
in fact happened to the photon trajectory, despite appearances to the contrary; this will also be 
discussed further in a future work. 

We will now consider massive particles, which travel along timelike geodesies. The only dif- 
ference in this case is that one must now impose the timelike geodesic normalisation constraint 
g^v^^^'*' = 1 = L/4. Let us again begin by considering a particle travelling in the j-direction. In 



this case, (41 ) still holds, but the normalisation constraint becomes 4/^ = A + R^y^. This yields 



At ? = 0, the right-hand side is still constant to first-order, mirroring the massless case, so massive 
particles motion in the j-direction is also essentially unaffected by the pancaking. 



Finally, we consider motion of a massive particle along the ^-direction. Again (44) holds, but 
the normalisation constraint is now 4/^ = 4 + R^x^, which yields 

^ = 7 = T . (48) 

^ Ri^\6Rl + K^ 
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Since Ri ^Ct near t = 0, for a constant C, on integrating we obtain 



x{t) 



C 



In 



(49) 



In this regime the denominator in the logarithm is approximately constant so massive particles will 
display qualitatively precisely the same behaviour as photons as they approach and leave the point 
of pancaking. 

Thus, in summary, we find that there are null and timelike geodesies that go smoothly through 
the Big Bang into a pre-Big Bang phase, but other such geodesies that spiral infinitely around a 
topologically closed spatial dimension. This issue does not seem to warrant too much concern for 
our model, owing to the periodic identification of these angular coordinates on the 3-spheres. This 
further supports our claim that this model is substantially better behaved than most conventional 
cosmological models. The issue of geodesic completeness will be further discussed in a future 
paper, but it is worth noting here that, even if the spacetime is incomplete, the singularity at the 
Big Bang can, at worst, be of quasiregular type [|43l as occurs also in Taub-NUT pni42| . 



V. REALISTIC BIAXIAL BIANCHI IX COSMOLOGY 



We now examine the viability of our model for providing a realistic cosmology. We begin by 
noting the argument of flU [51, that there is a natural, geometrical boundary condition on the uni- 
verse resulting from the need to match a Big Bang phase onto an asymptotic de Sitter phase within 
a particular type of conformal representation. Here a genuine cosmological constant A is being 
assumed, rather than a quintessence model. Using a conformal embedding and the symmetries of 
de Sitter space, a boundary condition is arrived at that the total elapsed conformal time should be 
equal to |. The condition is simply to demand that the future asymptotic de Sitter state be the fu- 
ture infinity surface of the conformal embedding. This amounts to imposing a boundary condition 
at temporal infinity, much like in (quantum) field theory. This singles out a particular flow line in 
{Q.,n,Q.A) space. 

In the context of a closed FRW model, it is shown in [|4l |5l that the conformal time constraint 
predicts within the correct range the degree of flatness of the universe, as well as the size of 
the cosmological constant. Moreover, the computed inflationary perturbations were shown to be 
consistent with WMAP data and could even account for the low-£ deficit. Extending this conformal 
time constraint to our Bianchi model, we find that this constraint can be satisfied by setting the 
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FIG. 1: Dynamics of the biaxial Bianchi IX model: evolution of the logarithm of the scale factors Ri and 
R2 in Planck lengths Ip versus log time (t in units of Planck time tp). 



values are for K set to 1, and are essentially a representative set, rather than having been fixed to 
get best agreement with current data.) These parameter values are all 'natural', but in order to fix 
the normalisation of the perturbation spectrum, the mass of the scalar field had to be rescaled and 
bo changes to a less natural value accordingly. This choice for the mass of the scalar field needs to 
be put in by hand for every model so does not constitute any unusual fine-tuning. 

This model looks like the universe that we observe in the sense that it comes out of a Big Bang- 
like state, followed by an inflationary phase and eventually reaches a state of steady expansion 
(see Figure [T]). There it can feasibly be matched onto a model of radiation domination followed 
by matter domination to recover our standard cosmology. Note, however, that this Bianchi model 
is interesting irrespective of our arriving at the particular parameter values above by using the 
conformal time condition. 

As mentioned earlier, the model exhibits no curvature singularity, with all physical quantities 
remaining finite through the Big Bang. The fact that the vanishing radius is odd in time results 
in parity inversion as we go through the Big Bang. The property that the energy density remains 
finite is contrary to the most common scenarios (although it was previously known that a massive 
scalar field can lead to a non-singular bounce in the special case of a closed FRW universe; see 



free parameters in our model to fc = 1, m 



~ 64000' 



ao = 1.2, bo = 18000 and fo = 13. (These 
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FIG. 2: Dynamics of the biaxial Bianchi IX model: evolution of the logarithm of the comoving Hubble ra- 
dius versus log time. Note that the curvature radius R{t) = (/?i/?2^3) ' > with corresponding Hubble function 
H{t) = ^ and associated Hubble radius H^. The parameters for the model are described in the text. During 
inflation, 1 / {RH) decreases with time, corresponding to accelerated expansion. Thus this model leads to a 
period of inflation lasting approximately for the period In? = 9 — 13.5. All quantities are in Planck units. 



14511 for more information on the history of such models). 

Figure [T] demonstrates that the scale factor i?i (t) approaches zero linearly as ? ^ 0, as expected 



from the series expansion in Section IV B Note that this behaviour results in a slope of 1 in the 
log-log plot as Int — oo. The other radii R2 = R3, as well as physical variables such as the 
scalar field and the scalar field energy density, tend to a constant at the Big Bang as shown in 
Figure [1} confirming our results from the series expansion. This shows that this model exhibits a 
pancake singularity. The relaxation of the assumption of isotropy makes a variety of singularities 
possible which generalise the pointlike FRW-singularity. For instance, pancake, barrel and cigar 
singularities are known to occur in Bianchi models HI 51 1371 . (In both barrel and cigar singularities, 
two scale factors tend to zero but the third increases without bound in cigar singularities, and ap- 
proaches a constant in barrel singularities.) It is worth noting, however, that a pancake singularity 
is atypical for Bianchi IX models. 

The evolution of the oblateness i?2(0 /Ri (0 (see Figure[T]), suggests that the evolution equations 
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favour the radii R2(t) = R3 (t) becoming similar to the radius Ri{t), in agreement with the need to 
be tolerably close to an FRW cosmology at late times. (This ties in nicely with the fact demon- 



strated in Section VI that the biaxiality of our model, i.e. i?2(0 = ^3 (0' ^1^° stable to sufficiently 
small perturbations of the type 7?3(?)/i?2(0 = 1 + thereby allowing an almost-FRW model 
in which all axes are similar.) Note, however, that Figure [T] demonstrates that the universe was 
significantly anisotropic until at least \nt ^ 10. The scalar field seems to isotropise the universe 
and subsequently drive inflation. (For related literature on isotropisation, see [|46ll47ll48ll49l[50l .') 

We can define an averaged scale factor by R{t) = {RiR2R3)^ , with which we can associate an 
averaged Hubble function H{t) = R/R, as usual, and thereby define a comoving Hubble radius 
1 / (RH). Figure |2] shows the evolution of the comoving Hubble radius 1 / (RH) for our model. 
Inflation is a period of accelerated expansion, which results in the comoving Hubble radius 1 / (RH) 
decreasing with time. From the plot we can infer that for our model a period of inflation does occur 
and lasts approximately from In? 9 to In? 13.5. During this inflationary period ln(i?) increases 
from about 10 to 55 (c.f. Figure [1]), corresponding to roughly 45 e-folds of inflation. 

In order to produce the fluctuation spectrum observed in our own universe, it is thought that 
the present Hubble horizon must have been inflated by at least about 40-50 e-folds, that is, 40-50 
e-folds between the time where the present Hubble scale exited the horizon and the end of inflation 
[|5Tll52ll . This minimal number of e-folds would lead to a universe that departs from spatial flatness 
by an amount that is just visible today ll53l . If more time was spent on the inflationary attractor, 
the universe would be driven closer to flatness and might in fact become indistinguishably close to 
spatially flat. 

Since our model produces around 45 e-folds of inflation, for agreement with observations, we 
would therefore require that the present horizon scale leave the Hubble horizon soon after the onset 
of inflation. This is in agreement with Uzan, Kirchner and Ellis' estimates [|52l and, in particular, 
also leads to a visibly closed geometry (cf. Starobinsky [|53ll ). 

Figures [T] and [2] show that the phases of isotropisation and inflation might have overlapped in 
the time period In? ^ 9 — 10 such that structure on the largest scales could have been laid down 
whilst the universe was still significantly oblate. In particular, we have just argued that the present 
horizon scale should have crossed the horizon around that time. So imprints of oblateness might 
actually be experimentally accessible to us. Such considerations, including perturbation analysis, 
CMB imprints, etc. will be described in future work. Clearly, a proper analysis must generalise 
standard inflationary cosmological computations [ISTl |54l to anisotropic cosmologies [|55ll . In 
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particular, there are subtleties regarding the correct choice of vacuum (generalisation of the Bunch- 
Davies vacuum), the definition of the canonical variables (generalisation of the Mukhanov-Sasaki 
variables), the discreteness of the eigenmodes of the Laplacian on the 3-sphere (as mentioned 
above, we use a continuum approximation) as well as the definition of k^, which arises since the 
direction of the wavevector kj starts to matter due to the anisotropy. Since inflation drives the 
spacetime closer to flatness, it might also be possible to approximate the spacetime by a Bianchi 
I model and quantise the perturbations in this approximation as an intermediate step. There has 
been recent progress in the cosmological perturbation theory of Bianchi I models, to which the 
interested reader may wish to refer mill El. Their results further show that the two gravitational 
wave polarisations do not necessarily have the same power spectra. The spectra also do not seem 
to reduce to the standard results in the limit of vanishing shear. At very early times, when the 
shear dominates over the curvature, the Bianchi IX model is also very close to a Bianchi I model 
- though not topologically, of course. Note also that since the issue of how geodesies propagate 
through the pancaking is as yet not completely resolved, we only consider perturbations 'this side' 
of the pancaking, rather than trying to track or match perturbations from both sides P31l56l[57ll . 

As a first approximation, however, one may consider the slight anisotropy during inflation as a 
perturbation to standard FRW results and neglect it to zeroth order. In this case, the definition of a 
comoving wavevector k becomes unambiguous. A detailed computation would ultimately have to 
validate this approach. In this approximation the power spectrum of the curvature perturbation is 
given by 

^M*) = 8.(|l)', (50) 
and the spectrum of the tensor perturbation by 

grav \ 



(51) 



where the right hand sides of these equations are to be evaluated when the corresponding comoving 
wavenumber k crosses the horizon. From this we can also define the tensor-to-scalar ratio 



grav 



(52) 

evaluated at some suitable low k. Note that there is a factor of 8;r premultiplying the standard 



result in equation (50). This is a consequence of the use of different conventions: G = 1 in the 



standard derivation, whereas here we have set SnG = K =1. 
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Non-flat universe models, of which Bianchi IX is an example, are more complicated than their 
flat counterparts in many respects. However, they have the decided advantage of having another 
length scale available at any time: their associated curvature scale R. As set out below, this allows 
one to compare length scales at different times directly, rather than having to use relations from 
the inferred history of the universe, such as the epochs of matter and radiation domination and 
reheating, the details of which are not well known. Note that current experiments do indeed allow 
for curvature contributions to Q. at the per cent level, making such computations consistent with 
observations. 

The comoving Hubble radius, and therefore the radius of the spatial sections, is related to the 
density parameter by 

The evolution of Q. for our model is shown in Figure [3] During inflation, is driven to unity 
from above, whereas after inflation Q. = \ becomes a separatrix rather than an attractor, and 
increases away from unity again. This provides a relationship linking time, the density parameter 
and the curvature radius (and thereby other length scales). We are particularly interested in when 
certain length scales cross the horizon. The first case of interest is when quantum fluctuations on 
different scales leave the horizon during inflation and thereby seed structure formation. Second, 
once inflation is completed, progressively larger scales then re-enter the horizon. The advantage 
of having the curvature scale available in non-flat geometries is that we can hnk the two times of 
horizon crossing for a particular length scale in a straightforward and exact manner. 

Consider some physical size do at the present time. (For wave modes in such a spherical 
universe, do will actually be quantised in units of lTtRo/{n^ — 1)^''^, where n is an integer. Here 
we will use a continuum approximation, however.) Scales grow commensurately with the scale 
factor during the expansion of the universe, such that 

^ = ^. (54) 

R Ro 

(Here, quantities without indices are evaluated at an arbitrary time, i.e. their dependence on time is 
left implicit, whereas the subscript denotes quantities evaluated at the present epoch.) Moreover, 
the size of the Hubble radius relative to the radius of the spatial sections is also changing over the 
course of cosmic history, as given by ([53]). Suppose do occupies some fraction xq of the current 



Hubble radius Hq ^ 



xo = — Y- (55) 

^0 
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FIG. 3: Evolution of Q.. During inflation the total density is quickly driven to 1. 



The ratio x — dH will of course in general also change over time and obeys 

R RH 

X = dH = dQ—H = xq—— . 



(56) 



Horizon crossing occurs when x=\, i.e. the scale with a size the present time left the horizon 
(i.e. Hubble length) at a time t given implicitly by 

RqHq ^^{t) - 1 



Xq 



(57) 



RH V^^o- 1 ' 

where is the present value of the density parameter. As a simple example, the scale re-entering 
the horizon at the moment is, of course, the current Hubble horizon itself (for which xq = 1). Thus, 



from (57 ), the present Hubble radius left the horizon at a time t when Q.{t) = Q.q during inflation. 
So for comparison with observations, we are interested in scales that left the horizon after this 
epoch. 

The WMAP 3-year resuks [6| quote 1.01 1 ±0.012. We shaU therefore assume that Hq ~ 1.01. 
By the argument above, the time t at which the present Hubble radius left the horizon is given 



by the solution of Q.{t) = 1.01 during inflation. This is, by virtue of equation (53 1, equivalent to 
— \n{RH) = — In 10 ~ —2.3, which, referring back to Figure|2} occurs at ln(?) ^ 10.4. Comparing 
with Figure [T] this corresponds to about one e-fold of inflation. This still leaves 44 e-foldings 
before the end of inflation, which is sufficient to provide a perturbation spectrum consistent with 
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FIG. 4: Evolution of H /(p during inflation. The function displayed is — 10 // /0 as a function ofln{t ftp). 
This function determines the amplitude of the curvature perturbation. 

experiments, in agreement with the estimates in [52 1. Intriguingly, however, we can also extract 
the level of oblateness at that time from Figure [T] as ~ 0.2%. This is to be regarded as significantly 
anisotropic, which means that structure on the largest scales would have been laid down when the 
universe was still oblate. This offers the exciting prospect of a possible experimental detection of 
imprints from such a time. 



Having clarified how to evaluate the quantities given by ( |50| ) and ( [51] ), we note that the power 
spectrum of the curvature perturbation is controlled by the quantity H^/ <j), plotted in Figure[4j The 
curvature perturbation spectrum can now be computed from the joint knowledge of the evolution of 
Q. (Figurepl) and //^/ (p (Figured). A physical wavenumber kp Q today is the inverse of the physical 



length scale do and from ([55]) we find the physical wavenumber today that left the horizon at the 

1 



time t from ( 57 ) is given by 



xo_ 

kp,o Hq ' 



or, in comoving terms. 



1 



^0 

RqHq 



(58) 



(59) 



Thus for a given comoving wavenumber k, equation (59) yields the corresponding value of xq 
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FIG. 5: The power spectrum of the curvature perturbation. 4n^l0^ ^_^{k) is plotted as a function of 
ln(l/jco) = ln(/c) + constant. The grey line is a fit to the spectrum with an exponential cutoff proposed 
by Efstathiou. The parameters of the model and the cutoff are described in the text. 



which is then substituted into ( [57] ) to obtain the time t when that scale left the horizon such that 

1 



v/ao-i = Va(0-i. 



(60) 



Figure [5] shows the curvature power spectrum as a function of \n{\/xo) = \n{k) + constant (= 
ln{kp o) + constant). We again stress that, firstly, due to the isotropisation we are using an effective 
FRW-description, such that the direction-dependence of the modes is negligible to zeroth order. 
Secondly, the eigenmodes in a closed universe are discrete - the spectrum shown here is, strictly 
speaking, a continuous approximation to the underlying discrete spectrum. It is essentially the 
direction-averaged power spectrum that is an enveloping function to the true quantised spectrum. 

The spectrum exhibits a sharp cutoff at low k, which is a consequence of the quantity H^/<l) 
turning over in the range t ^ 9.5 — 10.5 (which is not the case in conventional flat FRW models 
with straightforward power laws). Such an exponential cutoff of the form 



<^^(/t) = a{l+b\n{k)){\ -exp{c{\n{k)+d)) 



(61) 



has in fact been argued for on phenomenological grounds by Efstathiou [|58l . Indeed, our pre- 
dictions agree rather well with this proposed exponential cutoff, cf. the grey line in Figure [5] It 
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is intriguing that our model fits phenomenological predictions, and, in particular, this could ac- 
count for the observed dip at low I in the CMB power spectrum. From the power law part of 
the spectrum, we can extract the spectral index as = 0.975, which is also broadly in agreement 
with observations [[6l. The tensor spectrum is qualitatively very similar to the scalar spectrum, 
and yields a tensor- to- scalar ratio of approximately r ~ 0.2 in agreement with current constraints. 
Such agreement with data is encouraging. However, we must again stress that these computations 
are to zeroth order and more detailed computations in the anisotropic setting are needed. 

The process of averaging the radii to form R= {RiR^^ is effectively a map from a deformation 
of an FRW model back to an FRW model. To zeroth order, we could then apply all the standard 
machinery for FRW universes. This is certainly a very sensible way of achieving an FRW model 
from a Bianchi model, but is to some extent arbitrary. Looking at two other FRW models derived 
from our Bianchi model is instructive. One could define an FRW model simply by using each of 
Ri and R2 in turn. Of course, that will lead to problems with singularities at early times, but for 
the purpose of calculating fluctuation spectra one could use the numerical value of either radius 
in our Bianchi model just before inflation as an initial condition for an FRW model. It then turns 
out that in a model that has the appropriate value of i?i as a starting point for an FRW universe at 
the beginning of inflation, the quantity which controls the magnitude of curvature perturbations, 
H'^/<j), does not turn over. The curvature spectrum in this model therefore no longer exhibits an 
exponential cutoff of the above mentioned form. A model based on the other radius, R2, however, 
does have a turnover similar to the one in Figure |4] and therefore an exponential cutoff in the 
curvature perturbation spectrum. This then explains how the unusual feature of an exponential 
cutoff arises via the averaging over the two differently behaved radii. 

In summary, the model leads to isotropisation, necessary for compatibility with standard cos- 
mological models at late times, as well as inflation, which accounts for structure formation. The 
perturbation spectrum we predict meets current constraints on the spectral index and tensor-to- 
scalar ratio, and offers an explanation for the dip at low i in the CMB power spectrum. There is 
also the intriguing possibility of imprints of early oblateness on structure formation. We suggest 
that the universe could have been 0.2% oblate at the time when the present Hubble radius left the 
horizon. 
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FIG. 6: Dynamics of the full triaxial Bianchi IX model: The natural logarithm of the ratio of the nearly 
degenerate radii is plotted as a function of log time (in Planck units). 

VI. PERTURBATIONS AROUND THE AXISYMMETRIC CASE 



Given the remarkable properties of this axi symmetric model, it is important to know how stable 
it is to perturbations. If a small fractional perturbation R3{t) /R2{t) = 1 + 5{t) in the initial radii 
evolved to universes in which ^2(0 and Rsit) were vastly different, the axisymmetric case would 
hardly be a viable model for cosmology. There are several ways in which we can study the stability, 
both numerically and analytically. 

Setting Rsit) = R2{t){l + 5{t)) in the full triaxial equations yields the following dynamical 
equation for the fractional perturbation 5{t) 



[1 + 8)8 



4g4 16^3 {24Rl-4Rl) ^, 

-i^iH,+2H2)8 + ^^]8 
-{Hi+2H2)8. 



(62) 



One can either evolve this 5 (?) -equation, or alternatively one could solve the full triaxial equations 



(19)-(20) numerically, subject to the initial fractional perturbation 5(0) ^ 1. 

We chose the latter approach, i.e. numerical integration starting from a point close to the Big 
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Bang with a very small fractional perturbation 5(0) = 10^^ around the biaxial cosmology consid- 
ered in the previous section, keeping the other constants at their previous values. The numerical 
results show that as t ^ (ln{t) — > —oo) approaches zero, but the ratio of R3/R2 undergoes 
an infinite set of oscillations. (The evolution of the quantities corresponding to Figures [T] to |3] are 
indistinguishable from their biaxial counterparts and hence not shown again.) This is indeed what 
we would expect for the triaxial case (as cited previously, [|37i ; cf. Figure [6]). However, we are 
interested in what happens to these oscillations at large t rather than at early times. The results 
show that the oscillations cease and the ratio gets frozen in around the onset of inflation. This 
shows that for small perturbations a steady state is eventually reached where the ratio R3{t) /Riit) 
reaches a constant value. In particular, the resulting value of Rsit) /Riit) differs from unity by an 
amount that is similar to the initial fractional perturbation. Hence the small perturbation does not 
lead ultimately to vastly different radii. 

Note that here we mean by ? = a different time from in previous sections. Until now, t = 
was naturally defined to be the time where pancaking occurs. However, perturbing and therefore 
going to the full triaxial equations will in general lead to an oscillatory singularity (see above, 
[|37ll ). Instead, we now take f = to mean the point where we impose the boundary conditions. 
These conditions are exactly the same those applied previously at the Big Bang, except that now 
the boundary condition for the perturbed radius Rj, is Rj,{0) = 7?2(0)(1 + 5(0)) = Z7o(l + 5(0)). 

It is worth noting that, as illustrated in Figure[6} the ratios of the scale factors become frozen-in. 
Thus, even if the values of the scale factors in each direction differ slightly, the Hubble functions 
Hi{t) in each direction will be equal. This then raises the question as to whether such a situation 
would be observationally distinguishable from an FRW universe at all. It does not seem to be in the 
flat case, but the fact that closed and open universes have an absolute distance scale associated with 
them via the curvature scale seems to suggest that it would be detectable. These considerations 
will also be described in future work. 

In fact the behaviour of the small perturbations around the biaxial case is reminiscent of the 
evolution of the conventional curvature perturbations generated from quantum fluctuations con- 
sidered earlier. These oscillate on subhorizon scales but get frozen in on superhorizon scales. 
They are often regarded as mini-FRW-universes that are locally over- (or under-) dense and evolve 
in another FRW-background model. It seems obvious that these can in general be anisotropic, so 
that Bianchi models, in particular Bianchi IX, might be good for describing cosmological pertur- 
bations even when the background model is strictly Friedmann-Robertson-Walker. The different 
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behaviour on subhorizon and superhorizon scales is normally explained by saying that subhori- 
zon scales describe causally connected regions where 'there is time for local differences in matter 
distribution to affect the physics'. Superhorizon scales, however, are causally disconnected and 
frozen in. The cosmologically relevant scales that were previously subhorizon are stretched to 
superhorizon scales by inflation, and are now slowly coming back into the horizon. Therefore it is 
only natural to think of the perturbations in the Bianchi model in the same way - as stretched to 
superhorizon scales by inflation and frozen in. 

We can make this empirical connection slightly more quantitative. Consider the evolution of 
a curvature perturbation associated with a comoving wavenumber k in a flat FRW-universe [|54l . 
(Note that this is an excellent approximation for our model after the start of inflation.) This is 
given by 

'Ck+[§ + 2^+3Hy, + ^^^, = 0. (63) 

This equation is of harmonic oscillator form. However, amplitude and frequency are time- 
dependent. Thus it is reasonable that, in general, the motion will be oscillatory. At late times 
(after inflation), however, the effects of the scalar field driving inflation are negligible and the 
Hubble radius is large, so overall the damping term is negligible. Also the size of the universe is 
then large, so that the equation basically reduces to l^k = 0- Thus at late times the perturbation gets 
frozen in. 



Now consider perturbations around the biaxial case ( 62 1 in the late-time, isotropic limit Hi — H, 
Ri = R. This yields 

d + 3Hd + ^d=0. (64) 

Thus the two different kinds of perturbation actually obey very similar equations. They are prac- 
tically of the same form when the effect of the scalar field is negligible (after inflation). The fact 
that the terms involving the respective function and its first derivative are suppressed by the expan- 
sion of the universe explains why both types of oscillation become frozen-in eventually. From the 
equations, we would, however, expect the evolution of d to exhibit much less variation of ampli- 
tude and period with time, which is indeed what is observed numerically. Hence, at the level of the 
dynamical equations, it is reasonable that the two different sorts of perturbation look so similar. 
The deeper question of why they should obey such similar equations will be considered in a future 
work. 

One could, alternatively, argue that the more natural variable to consider is the fractional per- 
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turbation in the Hubble functions, rather than the radii. We choose to parameterise this fractional 
perturbation in an analogous way as H2,{t) = H2{t){l +h{t)). Its evolution is governed by the 
equation 

h = -^{Hi+H2)h^ 

1 / R^ 111 \ h 

2 V RjRj R\Rl R\Rl R\ Rj Rj J H2 
R] rI 1 1 \ 1 



R^R^ R\R'^ R^ R2 ) ^2 
From this equation it can be seen that, in fact, an expansion in which two axes have the same 
Hubble function is favoured and stable: for a small initial fractional perturbation the terms in 
powers of h are negligible, and the constant term vanishes when R2= R3, which implies that h = 



and, from (65 ), h = 0. This further supports the claim that this biaxial model is stable and hence a 



viable and interesting cosmological model. 



VII. BOUNCE SOLUTION 



In addition to the odd-parity, pancaking solution that we have discussed so far, there also exists 



an even-parity, bouncing solution to the biaxial evolution equations (25-28). In this solution, all 



radii are finite at the Big Bang and go smoothly through into a pre-Big Bang phase. The series 



solution ansatz akin to (30) is even under parity 

Ri (t) = ao + (^2t^ + ciAt^ + ■ ■ • 

R2{t)=R3{t)=bQ + b2t^ + Ht^ + ... . 

0(O=/o+/2?'+/4?' + ... 

(66) 

In contrast to the previous series solution, the Friedmann constraint is not automatically satisfied, 
and leads to an additional constraint determining uq in terms of Bq, fo, etc. Figure |7] shows the 
basic dynamics of this model. Both radii are constant and non-zero at the Big Bang, and the effect 
of the scalar field is again to inflate and isotropise the universe. For considerations concerning a 
presentation in terms of shear and curvature, see the appendix. This model might be interesting 
in its own right, and further details will be presented elsewhere. For previous work on bouncing 
solutions and their likelihood see, for example, 
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FIG. 7: Dynamics of the bouncing solution: evolution of the logarithm of the scale factors Ri and R2 in 
Planck lengths Ip versus log time (t in units of Planck time tp). 

VIII. CONCLUSIONS 

We have presented a novel scenario in which the effect of a scalar field and biaxial Bianchi IX 
geometry is to render the Big Bang much better behaved. Physical quantities remain finite at the 
Big Bang, and there is no curvature singularity at the time of pancaking. This is in contrast with 
singularities in flat and open FRW-models, and essential singularities in the full triaxial Bianchi 
IX case. In addition, the investigation of the behaviour of geodesies suggests that some physical 
observers can cross over and move smoothly into a pre-Big Bang phase. Some observers, however, 
may not be able to cross and this would appear to make the model incomplete, with the pancake 
being a quasiregular singularity (the type that also occurs in Taub-NUT). However, in the light 
of periodic identification of angular coordinates on the 3-sphere and the fact that the winding di- 
rection is the collapsing dimension, this might not be such a problem. The model also exhibits 
stabiUty under sufficiently small perturbations around biaxiahty. Though the model isotropises at 
late times, we show that structure on the largest scales could in fact stem from a time at which the 
universe was significantly anisotropic. Such behaviour is very desirable for a satisfactory cosmo- 
logical model, as a closed universe can be matched to the observed asymptotic de Sitter phase of 
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A- domination, whilst the scalar field can drive isotropisation, so as to account for the observed 
degree of isotropy of the universe, and inflation, thought to be needed in order to explain the ob- 
served perturbation spectrum. We predict a spectral index and a tensor-to-scalar ratio compatible 
with current constraints, as well as a dip at low multipoles of the CMB power spectrum, which is 
consistent with an exponential cutoff that has also been argued for on phenomenological grounds. 
We explain the apparent qualitative similarity between the two very different kinds of perturba- 
tions, namely curvature perturbations and perturbations around biaxiality, from their evolution 
equations. A separate, bouncing solution is also presented, and will be considered in more detail 
elsewhere. A more thorough analysis of many subtleties of the model, in particular the issue of 
geodesic completeness and a detailed comparison with Taub-NUT is in preparation. 
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APPENDIX: DYNAMICS IN THE 3 + 1 COVARIANT APPROACH 

Here we use the timelike eigenvector = [1,0,0,0] of the scalar field energy-momentum ten- 
sor as the 'velocity field' to perform a 3 + 1 split (see, for example, [61]) in order to present the 



equations ([25|)-(|28|) in a different parametrisation, which might be convenient for other applica- 
tions, and in particular separates more clearly the contributions of curvature and shear. Many 
quantities in the general 3 + 1 approach vanish for our particular case, such as the peculiar accel- 
eration, vorticity and anisotropic stress. The dynamical equations can therefore be recast simply 
in terms of the averaged scale factor, i.e. the volume expansion, 

R={RiR2R3)''^{RiRiy: (A.l) 

its associated Hubble factor, 

3H = Hi+H2 + H3^Hi + 2H2, (A.2) 
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and the shear scalar 



2 1 

^ =6 



{Hi-H2Y + {H2-H^Y + {H^-HiY ^^{Hi-H2)\ (A.3) 



where the expressions on the right denote the biaxial limit. The shear tensor a^v is defined as 
o^v = ^{fiiiy)^ where V is the fully orthogonally projected covariant derivative (for details, see 
ll6Tll ). The shear tensor has three spacelike eigenvectors characterised in a coordinate-free way 
by their eigenvalues. In the biaxial case, for instance, one is | {Hi —H2), whereas the other two 
are degenerate and have the value —4 {Hi —H2), whence the above shear scalar follows from 



a2 ^ icT^vCT^"- 



In this parametrisation, the Raychudhuri equation assumes the form 



3H + 3H^^-^{p+3p)-2o^+A. (A.4) 
The Friedmann equation can be reexpressed as 

where R^^^ is the Ricci scalar of the orthogonal 3-spaces, which here evaluates to 
r(^) = 2(^ + ^ + ^-^-^-^]^^(^- 

\ p2p2 p2p2 p2p2 p2 p2 p2 / p2 I p2 



(A.6) 



This uncovers the geometric significance of the corresponding term in the triaxial Friedmann equa- 
tion (21 ). We also find that the —H1H2 — H2H2 — //3//1 term just corresponds to 3H^ — (up to 
sign), and 3H^ + 2o^ = H^ +H2+ H^. This also correctly reduces to the standard FRW results 
in the isotropic limit, as do the equations themselves. The usual shear propagation for vanishing 
pecuhar acceleration and anisotropic stress is 

(7(^ V) + ^Ho^y = R^j^l - ^h^yR^^^ , (A.7) 

where R^^l and R^^^^ are the Ricci tensor and scalar of the orthogonal 3-spaces, respectively. How- 
ever, in the biaxial case, the evolution equation for the shear scalar suffices, and turns out to be 

a + 3Ha = -^i^^a^^ = _ \ _ ^^^^^ 



2g^' y3i?2 ^^2 



This can be seen upon contracting (A.7) with o^^ and noting that d^^y^a^^ = d^yO'^^ = 
2((7^v<7^^)' = 2ad, i.e. since one is contracting with something purely spatial, the projection 
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onto the 3-spaces prior to the contraction does not change the resuh. Furthermore, the term in- 
volving the Ricci scalar of the 3-surfaces just extracts the trace of the shear, which vanishes by 
definition. 

The equation of motion for the scalar field just retains its usual Klein-Gordon form 

m^^ + 3H<j) + ^=0. (A.9) 



Using the series solution (30) for the pancaking solution, we find that 

R{t) = {RiRj)'^ , (A. 10) 



to lowest order int. By the same token, the shear scalar behaves as 



as might be expected from the naive form of equation (A. 8 1 d + 3Ho = 0, which holds in the 
better-known flat and isotropic limits. Note that this does not contradict equation (A. 8 1. The 
lowest order term in the series expansion for a is a °^ Hi{t) °^ j which results in d — ^ to 
lowest order. This happens to cancel exactly the ^ term resulting from the product of c y and 



H{t) //i (?) oc such that the leading order of the combined o + 3Ho is actually the constant that 



is the leading term on the right hand side of (A. 8 1. Also note that the approximate time-dependence 



of the shear scalar does not involve any of the parameters that need to be set in our model, which 
is remarkable. The curvature of the orthogonal 3-space is constant across the pancake 

i?(3) = _^, (A.12) 

whereas the shear diverges, as shown. This means that in a pancaking Bianchi IX model the shear 
will always dominate over the curvature at early enough times, such that the model is actually well 
approximated by a Bianchi I model, though they are of course different topologically. 

Figure [8] depicts the motion of the pancaking model through the shear-curvature-parameter 
space. Note that since the shear diverges at the time of pancaking we have chosen to plot the 
model's trajectory in the (in a, R^^^^- space rather than the ^a,i?*^^)j -space itself (which will be 
appropriate for the bouncing solution). At early times, the shear diverges, whereas the curvature 
is constant, which accounts for the divergence in the top left of the plot corresponding to the time 
of pancaking. During inflation the curvature is quickly driven to zero and the universe isotropises 
i.e. the shear vanishes, which accounts for the late-time behaviour in the bottom right part. 
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FIG. 8: (a,/?^^^)-trajectory for the pancaking solution: evolution of the logarithm of the shear versus the 
Ricci scalar of the orthogonal 3-surfaces. 



Similarly, for the bouncing solution (66 1, we find that 



constant, (A. 13) 



and also the shear is constant (in fact vanishes) to first order 



2 4 /a2 ^2^.2 



<y =^[--7r] ^ (A.14) 

3 \ao bo J 

as might be expected, since all the dynamical variables have an extremum at the bounce. The Ricci 
scalar of the 3-surfaces is also constant 

i?(3)^ 4 (4-4V (A.15) 



bl Wo . 

This means that for a bouncing model, contrary to the pancaking case, the curvature will generi- 
cally dominate over the shear at the bounce. 

Figure |9] shows the motion of the bouncing model through the shear-curvature-parameter space. 
At early times, the shear vanishes, and likewise at late times. The shear therefore peaks at inter- 
mediate times in the bouncing model. The curvature shows the same behaviour as the pancaking 
solution in that it is constant at early times and the model is then driven to flatness during inflation. 
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FIG. 9: ((7,/?(^))-trajectory for the bouncing solution: evolution of the shear versus the Ricci scalar of the 
orthogonal 3-surfaces. 

So we can identify the bottom left hand area of the plot with the bounce, where the curvature is 
maximal and the shear vanishes. The shear then peaks at intermediate times whilst the curvature 
decreases. The bottom right hand corner again corresponds to the late-time evolution which is 
driven to isotropy and spatial flatness. Note that in both cases the numerical value of the initial 
curvature is exactly what we expect for the chosen parameter values. They appear to be nearly the 
same since for the particular models chosen here ao <^ &o- 



[1] T. S. Pereira, C. Piti-ou, and J.-P Uzan, JCAP 9, 6 (2007), 0707.0736. 

[2] C. Pitrou, T. S. Pereira, and J.-P Uzan, JCAP 4, 4 (2008), 0801.3596. 

[3] A. Giimriikfuoglu, C. Contaldi, and M. Peloso, JCAP 11, 005 (2007). 

[4] A. Lasenby and C. Doran, Phys. Rev. D 71, 063502 (2005). 

[5] A. Lasenby and C. Doran, AIP conf. proc. 736, 53 (2004). 

[6] D. Spergel et al., ApJS 170, 377408 (2007). 



35 



[7] A. V. Toporensky and V. O. Ustiansky, gr-qc/9907047 (1999). 
[8] S. Fay and T. Lehner, Gen. Rel. Grav. 37, 1097 (2005). 
[9] A. Coley and M. Goliath, Phys. Rev. D 62, 043526 (2000). 
[10] C. W. Misner, Phys. Rev. Lett. 22, 1071 (1969). 

[11] B. A. Belinsky, I. M. Khalatnikov, and E. M. Lifshitz, JETP 62, 1606 (1972). 

[12] B. A. Belinsky, I. M. Khalatnikov, and E. M. Lifshitz, JETP 60, 1969 (1971). 

[13] G. R R. ElUs and M. A. H. MacCallum, Commun. Math. Phys. 12, 108 (1969). 

[14] 1. D. Novikov and Y. B. Zel'dovich, Annu. Rev. Astro. Astrophys. 11, 387 (1973). 

[15] G. F. R. EUis, Gen. Relativ. Gravit. 38, 1003 (2006). 

[16] J. M. Heinzle and C. Uggla, ArXiv e-prints (2009), 0901.0776. 

[17] J. M. Heinzle and C. Uggla, ArXiv e-prints (2009), 0901.0806. 

[18] J. M. M. Senovilla, Phys. Rev. Lett. 64, 2219 (1990). 

[19] E. Ruiz and J. M. M. Senovilla, Phys. Rev. D 45, 1995 (1992). 

[20] F. J. Chinea, L. Fernandez- Jambrina, and J. M. M. Senovilla, Phys. Rev. D 45, 481 (1992). 

[21] L. Fernandez- Jambrina and L. M. Gonzalez-Romero, Phys. Rev. D 66, 024027 (2002). 

[22] D. S. Goldwirth and T. Piran, Physics Reports 214, 223 (1992). 

[23] S. Kanno, M. Kimura, J. Soda, and S. Yokoyama, JCAP 8, 34 (2008), 0806.2422. 

[24] B. Himmetoglu, C. R. Contaldi, and M. Peloso (2008), 0809.2779. 

[25] L. Bianchi, Mem. della Soc. It. delle Scienze (detta dei XL) (3) 11, 267 (1897). 

[26] A. H. Taub, Ann. Math. 53, 472 (1951). 

[27] O. Heckmann and E. Schiicking, Gravitation (Wiley, New York, 1962). 

[28] J. Wainwright and G. Ellis, Dynamical Systems in Cosmology (Cambridge University Press, 1997). 
[29] C. B. Collins and S. W. Hawking, The Astrophysical Journal 180, 317 (1973). 
[30] G. F R. EUis and M. A. H. MacCallum, Commun. Math. Phys. p. 1072 (1968). 
[31] S. W. Hawking and G. F. R. ElUs, The Large Scale Structure of Space-Time (Cambridge University 
Press, 1973). 

[32] H. Stephani et al.. Exact Solutions to Einstein's Field Equations, 2nd edition (Cambridge University 
Press, 2003). 

[33] C. Doran, A. Lasenby, and J. Lasenby, Uncertainty in Geometric Computations (Kluwer Academic 
PubUshers, Boston, 2002), chap. Conformal Geometry, EucUdean Space and Geometric Algebra, pp. 
41-58. 



36 



[34] C. Doran and A. N. Lasenby, Geometric Algebra for Physicists (Cambridge University Press, 2003). 

[35] J. D. Barrow, Nature 272, 211 (1978). 

[36] J. D. Barrow, Nucl. Phys. B 296, 697 (1988). 

[37] M. P. Ryan and L. C. Shepley, Homogeneous Relativistic Cosmologies (Princeton University Press, 
1975). 

[38] K. Becker, M. Becker, and J. H. Schwarz, String Theory and M-Theory (Cambridge University Press, 
2007). 

[39] T. Ortm, Gravity and Strings (Cambridge University Press, 2004). 

[40] C. V. Johnson, D-Branes (Cambridge University Press, 2003). 

[41] D. A. Konkowski, T. M. Helliwell, and L. C. Shepley, Phys. Rev. D 31, 1178 (1985). 

[42] D. A. Konkowski and T. M. Helliwell, Phys. Rev. D 31, 1 195 (1985). 

[43] G. F. R. ElUs and B. G. Schmidt, Gen. Rel. Grav. 8, 915 (1977). 

[44] D. A. Konkowski and T. M. Helliwell, in The Tenth Marcel Grossmann Meeting. On recent devel- 
opments in theoretical and experimental general relativity, gravitation and relativistic field theories, 
edited by M. Novello, S. Perez Bergliaffa, and R. Ruffini (2005), pp. 1829-+. 

[45] C. Gordon and N. Turok, Phys. Rev. D 67, 123508 (2003). 

[46] R. M. Wald, Phys. Rev. D 28, 2118 (1983). 

[47] J. D. Barrow, Phys. Lett. B 187, 12 (1987). 

[48] Y. Kitada and K. I. Maeda, Phys. Rev. D 45, 1416 (1992). 

[49] Y. Kitada and K. I. Maeda, Class. Quantum Grav. 10, 703 (1993). 

[50] P Salucci and R. Fabbri, Nuovo Cimento B Serie 77, 62 (1983). 

[51] A. R. Liddle and D. H. Lyth, Cosmological Inflation and Large Scale Structure (Cambridge University 
Press, 2000). 

[52] J.-P Uzan, U. Kirchner, and G. F. R. Ellis, MNRAS 344, L65 (2003). 

[53] A. A. Starobinsky, Cosmoparticle Physics (Frontiers, 1996), chap. Spectrum of initial perturbations in 

open and closed universes, p. 43. 
[54] M. P. Hobson, G. P. Efstathiou, and A. N. Lasenby, General Relativity: An Introduction for Physicists 

(Cambridge University Press, 2006). 
[55] K. Dimopoulos, D. H. Lyth, and Y. Rodriguez (2008), 0809.1055. 
[56] S. Gratton, J. Khoury, P J. Steinhardt, and N. Turok, Phys. Rev. D 69, 103505 (2004). 
[57] J. K. Erickson, S. Gratton, P J. Steinhardt, and N. Turok, Phys. Rev. D 75, 123507 (2007). 

37 



[58] G. Efstathiou, MNRAS 343, L95 (2003). 

[59] J. D. Barrow and R. A. Matzner, Phys. Rev. D 21, 336 (1980). 

[60] F. T. Falciano, M. LiUey, and R Peter, Phys. Rev. D 77, 083513 (2008), 0802.1 196. 

[61] G. F. R. Ellis and H. van Elst, Cosmological models, Cargse Lectures (1998). 



38 



